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Abstract. We study the evolution of local event- by-event deviations from smooth 
average fluid dynamic fields, as they can arise in heavy ion collisions from the 
propagation of fluctuating initial conditions. Local fluctuations around Bjorken flow 
are found to be governed by non-linear equations whose solutions can be characterized 
qualitatively in terms of Reynolds numbers. Perturbations at different rapidities 
decouple quickly, and satisfy (after suitable coordinate transformations) an effectively 
two-dimensional Navier-Stokes equation of non-relativistic form. We discuss the 
conditions under which non-linearities in these equations cannot be neglected and 
turbulent behavior is expected to set in. 

In recent years, hadronic transverse momentum spectra and their azimuthal 
dependence with respect to the orientation of the reaction plane have provided tight 
constraints on the fluid dynamic model of ultra-relativistic heavy ion collisions. With 
the first data from the LHC, and with refined analyses of RHIC data, higher order 
flow coefficients v$, v±, v$ and Vq are now starting to complement the measurements 
of elliptic flow t>2- This provides access to qualitatively novel features of the collision 
dynamics. In particular, since event-averaged initial conditions of heavy ion collisions are 
by construction symmetric with respect to the reaction plane at mid rapidity, the recent 
measurements of non- vanishing odd harmonic coefficients V3, v& provide unambiguous 
evidence for the relevance of event-by-event fluctuations in the fluid dynamic evolution. 
Remarkably, at least some of the models currently used to specify initial conditions of 
fluid dynamic simulations can provide naturally for initial event-by-event fluctuations of 
the phenomenologically required size [1] . First studies of the fluid dynamic propagation 
of such geometric initial state fluctuations have resulted in marked improvements in the 
comparison of fluid dynamic simulations with data [2]. 

Motivated by these developments, we are studying here the evolution of local event- 
by-event deviations from smooth average fluid dynamic fields within the expanding 
geometry characteristic for ultra-relativistic heavy ion collisions. We are for instance 
interested in how a 'primordial' spectrum of such fluctuations evolves, and which 
modes of the fluctuation spectrum are damped on which time scales. Given that fluid 
dynamic fluctuations measure deviations from equilibrium, we hope to gain in this 
way novel access to the basic problem of how equilibration can proceed efficiently in 
heavy collisions. We also wonder to what extent it is justified to propagate primordial 
fluctuations in a linearized ansatz that by the nature of its approximation leaves no 
room for the development of turbulent phenomena. Such an approach is well-motivated 
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in the treatment of primordial fluctuations in cosmology, but as we shall discuss in the 
following, the scales and flows in heavy ion collisions can support a qualitatively different 
conclusion. 

To be specific, we consider local fluctuations in the fluid dynamic velocity 5u^ and 
energy density Se on top of average fields u M and e that satisfy Bjorken's scaling solution. 
We work in light cone coordinates, where m m = (1,0,0,0) and e = eBj(r ) (t /t) 4//3 . In 
general, one can split at fixed time r an arbitrary velocity field m m = u M + 8u^ into an 
irrotational part represented by the divergence 

= d lU l + d 2 u 2 + d y u v , (1) 

and a solenoidal part represented by the vorticity field 

oji = rd 2 u y d y u 2 , u>2 = —dyU 1 — rdiu v , u 3 = diu 2 — c^w 1 . 

T T 

Because of the normalization u^u^ = — 1, there are only three independent components 
Svj which we choose to span the transverse plane and rapidity, j = l,2,y. The energy 
density e = e + Se can be characterized in terms of temperature deviations from its 
Bjorken value 

d = ln(r/T Bj (r)). (2) 

Within this set-up, we have derived evolution equations for Su^ and Se from the 
relativistic viscous fluid dynamic equations for and e. 

We first comment on the linearized evolution equations for Su^ , Se. In this case, 
the equations for $ and d are coupled and describe essentially the propagation of sound 
(modulo modifications due to the expanding background, see Ref. [1]). The vorticity 
modes satisfy a diffusion-type equation of motion that can be solved directly. In Fourrier 
space, denoting by k y the wave number conjugate to rapidity, 

[T\i -^s (^+fc 2 2 )(- 4/3 -o /3 )+^73^(-- 2/3 -o" 2/3 ) 
u j (T,k 1 ,k 2 ,k y ) = u j (r ,k 1 ,k 2 ,ky) i— j e 4t o 3r o 

(3) 

For this solution, we assumed a r- independent ratio rj/s that leads to a kinematic 
viscosity vq = rj/(e + p) = 7]/(sT) evaluated at time To. While the transverse vorticity 
components (h\ — hi — —2) fall off like l/r 2//3 for small wave- vectors or small kinematic 
viscosity u , the vorticity mode grows algebraically in this region (h s = 1). Therefore, 
at least in a linearized description, there can be initial fluctuations that do not attenuate 
within the phenomenologically relevant time scales of heavy ion collisions, see Fig. [TJ 

In general, a linearized formalism applies if perturbations are small 5e/e 1, d 1 
and if Reynolds numbers are small, 

Re = u T I (e + p) = u T lsT 

Here ut is a characteristic velocity of fluctuations in the transverse direction and I 
is the typical length scale over which it changes significantly. (For motion in the 
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Figure 1. Vorticity amplitudes ai3(r)/w3(lfm/c) (solid lines) and wi(r)/wi(lfm/c) 
(dashed lines) for wave vectors with k2 — k y — 0, and k\ — 5fm _1 , k\ = lOfm -1 and 
k% = 30 fm" 1 , respectively. Viscosity has been chosen as v{t)/t = 10 -3 at r = lfm/c. 




rapidity direction, there is a different Reynolds number [4].) If the Reynolds number 
is large, then turbulent flow sets in and a linearized formalism fails. In this case, the 
smallness of perturbations around an expanding background still leads to the important 
simplification that the fluid can be treated as compression-less, $ = 0. This condition 
does not mean that there are no sound waves present, but that the coupling between 
sound and turbulence becomes negligible. The formal criterion is a small Mach number 

\fuiU 1 + U2U 2 + u v u y , , 

Ma = y — < 1 , (5) 

where cs is the velocity of sound. It is conceivable that the condition Ma <C 1 for 
compression-less turbulent flow is realized in heavy ion collisions. 

In terms of a rescaled time t = 3r 4 / 3 /(4r 1 ^ 3 ), temperature field d = (r /r) 2//3 (i and 
velocity Vj = (ro/r) 1 / 3 ?/^, we find that the non-linear fluid equations (j = 1,2, y) take 
the form 

fyvj + ^2 v md m Vj + ^VydyVj + djd - z/ ( df + d\ + -^df \ vj = (6) 

m=l ^ ' 

with the solenoidal constraint d^Vi + c?2f2 + ^dyVy = 0. For late times, this becomes 
effectively a two-dimensional Navier-Stokes equation of non-relativistic form for a two- 
dimensional compression-less fluid. Therefore, the development of turbulent flow in 
heavy ion collisions can be discussed on the basis of an equation about which much 
is known already. In particular, for the three-dimensional Navier-Stokes equation, one 
knows that kinetic energy cascades in the case of turbulence from large structures in 
space to finer and finer ones where it is eventually dissipated. In contrast, in two 
dimensions, turbulent kinetic energy is subject to an inverse cascade from microscopic 
to more and more macroscopic structures. This plays a role for turbulent phenomena 
in the essentially two-dimensional layer of the Earth's atmosphere and it prompts us 
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Figure 2. Illustration of turbulent kinetic energy as a function of the wave-number 
for freely decaying turbulence according to Batchelors scaling theory. 



to wonder whether ([6]) allows for a scale amplification mechanism of some fluctuating 
modes in heavy ion collisions. 

One way to characterize the inverse cascade is to study how the kinetic energy in 
a fluid, 



is distributed over wave vectors E(k) as a function of time. Based on the theory for 
turbulence developed by Kolmogorov in three and Kraichnan in two dimensions, there is 
in particular a scaling theory for freely decaying two-dimensional turbulence that leads 
to E(t,k) = A 3 1 h(k Xt), where h(x) is conjectured to be universal [3J. In this case, 
kinetic energy accumulates at small wave- vectors at late times, see Fig. [2j 

The scaling theories of Kolmogorov, Kraichnan and Batchelor address turbulence 
at very large Reynolds number. Estimates of Reynolds numbers realized in heavy ion 
collisions are usually of the order Re ~ s/i] = O(10). This is too small for applying 
results from fully developed turbulence, but it is sufficiently large to motivate a search 
for the onset of turbulent phenomena involving fundamental quantum fields in heavy 
ion collisions. This is a report on work in progress. We are currently exploring possible 
signatures for the onset of turbulent phenomena in measurements of single inclusive 
hadron spectra and two-particle correlation functions. 
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